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We study the intrinsic orbital magnetization (OM) in antiferromagnets on the distorted face- 
centered-cubic lattice. The combined lattice distortion and spin frustration induce nontrivial k- 
space Chern invariant, which turns to result in profound effects on the OM properties. We derive 
a specific relation between the OM and the Hall conductivity, according to which it is found that 
the intrinsic OM vanishes when the electron chemical potential lies in the Mott gap. The distinct 
behavior of the intrinsic OM in the metallic and insulating regions is shown. The Berry phase 
effects on the thermoelectric transport is also discussed. 
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I. INTRODUCTION 



The orbital magnetism of Bloch electrons has been an outstanding problem in solid state 
physics, and attracted renewed interest due to the recent recognition that the Berry 

phase effect plays a very important role in the orbital magnetism. The issue was carried out 
in the powerful semiclassical formalism , in which the Bloch electron for n-th band is 
treated as a wave packet \w n (r c , k c )) with its center (r c , k c ) in the phase space. The orbital 
magnetic moment characterizes the rotation of the wave packet around its centroid and is 
given by m n (k c )=^^-(w ra |(f — r c ) x v\w n ), where (— e) is the charge of the electron and v is 
the velocity operator. By writing the wave packet in terms of the Bloch state, one obtains 
(k c is abbreviated as k) 

m n (k) = -z(e/2fr)(V k u nk | x [# k - eJ2]|Vk«nk>, (1) 

where |w n k) is the periodic part of the Bloch state with band energy e^, and is the 
crystal Hamiltonian acting on |« n k). Equation ([1]) can be alternatively derived by taking 
the differential of the electron energy, which within first order in the perturbative magnetic 
field B turns to be e n k=e^l — m n (k)-B, with respect to B. It was further found 1] that the 
presence of a weak magnetic field B will result in a modification of the density of states in 
the semiclassical phase space, d 3 k — > <i 3 k(l + eB-fi n /fr), where O n (k) = i(VkUnk| x |Vk^ n k) 
is the Berry curvature in fc-space. Due to this weak-field modification, a quantum-state 
summation J2k 0(k) of some physical quantity C(k) should be converted to an integral 
according to J d 3 \t(l + eT$-Q n /h)0(\t). Additional thermodynamic average over Bloch bands 
should be included at finite temperature. Therefore, the total free energy for an equilibrium 
ensemble of electrons in the weak field may be written as 

F = ~ I rf3k i 1 + Ji B " n «( k )) + e^-^]. (2) 

n 

where /i is the electron chemical potential and j3 = 1/ksT. The equilibrium orbital mag- 
netization (OM) density is given by the field derivative at fixed temperature and chemical 
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potential, M = — (<9F/<9B)^ T , with the result 

M = J2 I d 3 km n (k)/ n (3) 
+ \ E / rf3k ^«(k) In [1 + e^~^] 

n 

= M C + M n , 

where f n is the local equilibrium Fermi function for n-th band. In addition to the conven- 
tional term M c in terms of the orbital magnetic moment m n (k), the extra term in Eq. 
([3]) is a Berry phase effect and exposes a new topological ingredient to the orbital magnetism. 
Interestingly, it is this Berry phase correction that eventually enters the thermal transport 
current [3j. At zero temperature and magnetic field the general expression is reduced to 

/Mo , p x 

d 3 k (m n (k) + -fi n (k) i(i ~ £nk]J , (4) 
n 

where the upper limit means that the integral is over states with energies below the zero- 
temperature chemical potential (Fermi energy) fi . 

The Berry phase effect on orbital magnetism was until now partially presented by very 
few studies. Recent observation of the anomalous Nernst effect (ANE) in CuCroSeA-^Br^. 

n n 

compound |6j was attributed [3J to the manifestation of the Berry phase effect in the OM. 
Also the orbital magnetism was recently studied by use of two-dimensional (2D) Haldane 
model and ferromagnetic kagome lattice with spin chirality Q, I^L These two models are 
rare examples to show the zero-field quantum Hall effect (QHE) 9|, ]l0|. From Ref. 8] one 
learns that the Berry phase effect causes the OM to display different behavior in metallic 
and insulating regions. This difference may be explained in parallel with Haldane's recent 
finding of the Berry phase effect in the intrinsic Hall conductivity [including QHE and 
anomalous Hall effect (AHE)]. 

The objective of the present paper is dual. First we remark that the Berry phase effect 



on the orbital magnetism has been included in the well-known Kubo-Streda formula [12 ] . 
Therefore, a full quantum-mechanical linear response theory of the OM can be developed as a 
useful complement of the semiclassical formalism, although the latter looks more elegant and 
practical for calculation on clean samples. Then we examine the 3D problem by studying the 
orbital magnetism in antiferromagnets on the distorted face-centered-cubic (fee) lattice. The 
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results reveal that a general "topological orbital magnetism theory" that takes into account 
Berry phase effect must now be developed. 

The paper is organized as follows. In the next section, we address that the intrinsic OM 
given in the semiclassical formalism is consistent with the well-known quantum-mechanical 
Kubo-Streda formula in the clean-sample limit. Section III describes the physical model 
that is used in this work. The topological property and the consequent intrinsic Hall effect 
associated with the model are also given in this section. In Sec. IV, we present a detailed 
study of the properties of the OM and its effects on transport response in antiferromagnets 
on the distorted fee lattice. Finally, in Sec. V we present our conclusions. 



II. KUBO-STREDA FORMULA OF THE ORBITAL MAGNETIZATION 

Due to the above mentioned modification of the density of states, the particle number in 
the weak magnetic field (say, along z-axis) is given by 

N{B, n) = Y,f d " k i 1 + /»• ( 5 ) 

n 

It is easy to see that a link between Eq. (j3J) and Eq. (jSJ) is y^jf'J = (f^)^' which is 
nothing but the usual thermodynamic Maxwell relation and therefore should be free from 
the weak-field limit used by the semiclassical approach. Thus the zero-field OM is given by 
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On the other side, the integrand in Eq. can be written in terms of Kubo-Streda 
formula for electrons as follows 

_ l _ dN(B, /i) 

0~xy | n — 0~ xy | fj, e , {( ) 



where a xy \^ is the Hall conductivity and 

<U = ^ / de^^-Tr[v x G + (e)v y 5(e - H) - v x 5(e - H)v y G-{e)\. (8) 

Here G ± (/i) = Hm J? _ >0 +(A t— Hiir/)' 1 is the operator Green function and v a is the velocity op- 
erator. In Bloch-state representation the trace in Eq. (jSJ) is equivalent to Ylmk^\ u nk\i.' ' ' )\ u nk) 
with the Hamiltonian transformed to ifk = e t]l ' r He~ l]i ' r and the velocity to v a (k) = -^^ K . 
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Replacing in Eq. © by -^f 1 and using the completeness relation of the Bloch 

states, J2nk l w nk)(M„,k|=l, one has 

-,2; 



I a L\^'d^' = -i^TT hm \ de/(e,/i) x V{— — ^7^-(Mnk|^(k)|M n 'k')( M n'k'|^(k)|M nk ) 
8{e- £ nk ) 



n,k 
n',k' 



■ («nk|%(k) | M„'k' ) («n'k' \V X (k) | U nk ) } 



e - £ n 'k' - ^ 

e 2 v Im { ^ life I u «' k ' ) ^™' k ' I Uf I u « k > ) 

TT> /(^nk,^) 

ft , £ nk — Sn'k' 



/(e^A*) • (9) 

P. . 1 _ /!_/ 

n,k 
n',k' 

By use of the identity 

~dK M = ^ Kk) + {£nk ~ Fk)l W (10) 

and after a transformation of k-sum to an integral, Eq. (Q is ready to be simplified as 

A comparison of Eq. (TTTT) with Eq. ([I]) immediately gives the following relation 

al^dn' = eJ2 f d 3 km z n (k)f n (12) 

= eMi 2) 
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where the second line is obtained by using the definition of M c in Eq. ([3]). Thus one finds that 
the semiclassical expression for M c is equivalent to the quantum-mechanical expression for 
(1/e) J M (Jxylu'd/j, 1 . Note that although a 1 is a Fermi-surface term, the quantity f M a^A^d^ 
is a Fermi-sea term and all the Bloch states below \x should be accounted when calculating 



the OM. On the other side, in the clean limit, the Kubo 



ormula for the Hall conductivity 



Cxyl/j, can be written in terms of the Berry curvatures Q z n 3], 

V X y\» = -jYl j d * k fnK- (13) 

n 

From Eq. ffl3l) one has 

f °xy\»'dn' = j ^k^n(k)ln[l + e^-^], (14) 
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which is exactly the semiclassical expression for the Berry phase term eM^ in Eq. ([3]). 
Thus the semiclassical OM in Eq. (J3J) can actually be written as a Kubo-Streda formula: 




The other two components M. x and Ai y are given in a similar manner. This equivalence 
between the semiclassical and quantum-mechanical description for the OM is only valid in 
the intrinsic region and will break down when the impurity scattering effect is included. 
Thus while the semiclassical formula of the OM is more suitably employed to study the 
intrinsic property of the OM, the Kubo-Streda formula must be used when one takes into 
account the impurity scattering. Another aspect is that the Kubo-Streda formula is valid in 
arbitrary strength of the external magnetic field, while the semiclassical formula only works 
in the weak field limit. In some special cases, for example, when one wants to know the edge 
state effect on the OM in a finite-size sample in a strong magnetic field B HQ , a more 
apparent approach can be employed by directly calculating the total free energy and then a 
finite-field OM is obtained by a 5-derivative of the free energy. 



III. SPECIFIC MODEL AND 2D CHERN NUMBER 



Now we focus our attention to the properties of the OM in a specific spin-frustrated 
system. As in Ref. [16j, the model we used describes the chiral spin state in the ordered 
antiferromagnet (AF) on the three-dimensional fee lattice. The AF on the fee lattice is a 
typical frustrated system, and nontrivial triple-Q_spin structure with finite spin chirality 
has been revealed by band structure calculation 17] and observed in experiments Q, [l9| . 
The anomalous behaviors in the fee AF were also observed. For example, there occurs 



mysterious weak ferromagnetism in N1S2 below the second AF transition temperature 20]. 
The Hall conductivity in this material is also large and strongly temperature dependent 21]. 
In Co(S a; Sei_. :r )2, the AHE is enhanced in the intermediate x region, where the nontrivial 
magnetism is realized jzj. 

The triple-Q spin structure on face lattice is shown in Fig. 1. Here the lattice points 
are divided into four sublattices with different local spins S a (a = 1,2,3,4) on them. The 
AF nature requires ^2 a S a =0. The minimization of the 2-spin exchange interaction energy 
cannot uniquely determine the sublattice spin orientation. The inclusion of higher order 
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FIG. 1: (a) Triple-Q spin structure on fee lattice, (b) Relation of the 4-spin moments S a (a = 
1,2,3,4). 

(4-spin exchange) interaction if 4 = J 4 ^2 a ^ b {S a -Sb) 2 with positive J4 gives the ground-state 
spin configuration 23, M as £l=(^, 7=, 4=), &=(^, ^=(-7%, 73, -75), and 



^4 = ( — 75 > — 73' 73); wnere each direction corresponds to the four corners from the center 
of a tetrahedron [Fig. 1(b)]. The effective Hamiltonian for the hopping electrons strongly 
coupled to the mean-field effective magnetic field caused by these local spins is given by 
H = J^NN^ij^l^j w ith ttj^ — t{Xi\Xj) = te mij cos^ 2 -. Here the spin wave function \xi) is 
explicitly given by \xi) = [cos y, e 1 ^ sin y] T , where the polar coordinates are pinned by the 
local spins, i.e., (Xi\Si\Xi) — \ (sin^j cos^j, sin^j sin0j, cos^j). $jj is the angle between 
the two spins Si and Sj. The phase factor can be regarded as the gauge vector potential 
aJr), and the corresponding gauge flux is related to scalar spin chirality Xijk=Si-(SjXS k ) 



251 ] . In periodic crystal lattices, the non- vanishing of the gauge flux relies on the multiband 
structure with each band being characterized by a Chern number. The Chern number 
appears as a result of the spin-orbit interaction and/or spin chirality in ferromagnets. In 
ferromagnets the time- reversal broken symmetry is manifest, while in AF the time-reversal 
operation combined with the translation operation often constitutes the unbroken symmetry. 
In the latter case, the nonzero Hall conductivity o xy is forbidden. However, when there are 
more than two sublattices and the spin structure is noncollinear, this combined symmetry 
would be absent and finite a xy is not forbidden 

The net spin chirality for the ideal fee AF lattice in Fig. 1 is zero, because the spin 
chiralities are the vector quantities and the sum of these four vectors is zero. However, 
when the lattice is distorted along the [1,1,1] direction, then the non-zero net spin chirality 
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^ X M r R 

71(0,1,0) 7t(1,1,0) Jl(1,1,1) 



FIG. 2: Energy spectrum along high-symmetry lines in the first BZ for (a) d=0, (b) d=0.2, and 
(c) d=-0.2. 

occurs. Following Ref. we express the distortion along the [1,1,1] direction by putting 
the transfer integral within the (1,1,1) plane as ti n tra=l, while that between the planes as 
ti nter =l—d. As the unit cell is cubic shown in Fig. 1, the first Brillouin zone (BZ) is cubic: 
[— -, ^] 3 . From now on, we set a=l. Then the Hamiltonian matrix for each k is given 
by 

/ e -*i/ a e*i/i / 3 \ 
eMr/ 2 e -i i/ 3 



7i 



-i— t 



(16) 



A e*Tr/ 2 / 



where f x =2{l - d) cos(^+^)+2 cos(-^+^), / 2 =2(1 - d) cos(^+^)+2 cos(-^+^), 
^3=2(1 — d) cos(^-+ %)+2 cos(— ^-+^)- In this Hamiltonian, the two lower bands are fully 
degenerate, £ik=£2k(k)=— yfi + /| + /|, while the two upper bands are also degenerate 



with £3 ) 4k=A//i + 7| + fi ■ The band structure along high-symmetry lines in the first BZ 
is shown in Fig. 2. At d=0 [Fig. 2(a)], the upper and lower dispersions touch along the 
edge of the BZ and comprises an assembly of the massless Dirac fermions (Weyl fermions) 
in (2 + 1)D. For d>0 (elongation along [1,1,1] direction), all the Weyl fermions along the 
edge open a gap and turn into the massive Dirac fermions [Fig. 2(b)]. Therefore the gap 
fully opens in the density of states centered at zero energy. For d<0 (suppression along 
[1,1,1] direction), all the (2 + 1)D Weyl fermions along the edges open the gap as in the 
case of d>0. However, there occurs two additional (3 + 1)D Weyl fermions [Fig. 2(c)] at 
(k x , k y , k z ) = ±D(1, 1, 1), where D = arccos (tzt) and ± correspond to the right- and left- 
handed chirality 26J. Thus unlike d>0, the gap does not fully open in the case of <i<0 due 
to the presence of a new single contact between the upper and lower bands in the BZ. The 
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normalized eigenvectors are given by 



1 



/i£ik+»/2/3 r -i§ /2£ik-i/i/3 n 



= lk 

/2 



Elk 



(17) 



|«2k) = — 7= ( e*i^-, e~*f 1, ze*i^- 

' ' -/2 V £ 2k' £2k' ' £2k 



|«3k) 



V2 



/l£3k+»/2/3 

^kv^Tif' 



— e 



E3k 



0. 



E3k 



|«4k) = (e*£^-, e _i t^3-, 1, i^t,^- 



£4k ' 



E4k ' 



t"4k 



The Berry curvatures for these four Bloch bands are derived to have the form 



Zfc nk 



with 



dkp dk 7 dk/3 dk 7 dkp dk 7 



(18) 



(19) 



where (a,/3, 7) represent a cyclic permutation of (x, y,z). 

Let us see the Hall conductivity of this system [16J with d 7^ 0. In the integer filling case, 
the zero-temperature Hall conductivity is a sum of Chern invariant 
bands, 

occu „ 

<J X y = (e 2 /h) V / —±C n (h s ), 
C n (k z ) is given by 



271 ] over occupied Bloch 
(20) 



where the 2D Chern number 

C n \kz) 



28 



1 

'27T 



dk x dk y Q z n (k) 



(21) 



-— / z ■ (V k x A„(k)) . 

Jf-vriTrl 2 



Here A n (k)=i(u r j c | VkU n k) is the Berry phase connection (vector potential) for n-th band. To 
proceed one may first transform the integral of Vk x A n over the first BZ to the line integral 
of A n along the BZ boundary by use of Stokes' theorem, and then apply the complex contour 
integration technique and residue theorem to sinusoidal functions. After a straightforward 
derivation, one obtains the non-zero Chern number, Ci{k z ) = —sga(g(k z )), C 3 —sga(g(k z )), 
where 

g(k 2 ) =2 + 2(1 - d) cos(k z + 2k P ) (22) 



and kp = arctan 



(d— 1) cos k z — 1 
(d— 1) sinfc z 



. It is easy to verify that for d>0, the value of g(k z ) is always 



positive, independent of k z . For d<0, the sign of g{k z ) depends on k z in such a way that 



g(k z )<0 for k z e (-D,D) and g(k z )>0 for k z G [-7T, —D) U (-D,7r]. Note that the present 
choice of the other two Bloch states |«2k) and |^4k) makes them to have no contribution to 
the Chern number. 

However, the above purely mathematical calculation of 2D Chern number is not favored 
by theoretical physicists, who would like to resort to the physical connotation that the vector 
potential A„ and gauge flux fl n are endowed with. Correspondingly, here we present this 
gauge-field analysis of the lower band £ik as an example. The value of 2D Chern number 
Ci(k z ), which is confined to the (k x ,k y ) subspace at fixed k z , is invariant under gauge 
transformation |M / lk )=e* </3l( ' k - ) |Mik), A'^k) = Ai(k) — Vk<£i(k), where <^i(k) is an arbitrary 
smooth function of k. If the gauge choice for |itik) is well-defined everywhere in the whole 
(k x ,k y ) subspace in the first BZ, then its Chern number C\(k z ) will obviously be zero. 
However, at point ko = (k z + 2kp, k z , k z ), one can find that the wave function |wik) in Eq. 
(I17p is ill-defined since both its denominator and numerator are zero at this point. This 
means that the used gauge cannot apply to the whole BZ and one needs to render a gauge 
transformation to avoid the singularity at ko. For this one transforms the wave function to 

\u') = -L( -W-H/i h -ilhe^-ihh o _ VIM ^ _ (23) 

The new eigenvector |w' lk ) recovers the well-defined behavior at k ; the new singularity 
brought about is at k = (k z ,k z + 2k P ,k z ). Thus according to the two different gauge 
choices, the BZ cross section at fixed k z is now divided into two regions V and V as shown 
in Fig. 3 (k z =0). The wave functions |wik) are used onto the region V, while \u' lk ) apply to 
V. Note that there remains some freedom in the division of the BZ. Because | it ik) and |w' lk ) 
are ill-defined only at k and k' , respectively, we are free to deform this division as long 
as k (k )^V(V'). This corresponds to the gauge degree of freedom 3, 30]. At keVflV, 
the two choices of wave functions are different by a phase factor |M / lk )=e i</3l ^ k - ) |wik), i-e., 
A' 1 (k)=Ai(k) - V k y?i(k), where 



jtpiQz) _ / fi + fi ^£ik + ihfi / 9zL x 

V fi + fih^ + ihfs 1 ] 

Thus one obtains the value of nonzero 2D Chern number for lower band £i(k) as follows 

C x [k z ) = ~(f [Ai(k) - A;(k)] -dk (25) 
= -— <b d^i(k) = -sgn^fc*)), 

^ JdV 
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FIG. 3: (Color online). Division of cross section (k z =0) of the first BZ into two regions V (red 
area) and V' (blue area). Note that the projection of the singularities ko and kg onto (k x ,k y ) 
subspace vary with k z . 

which is consistent with the explicit calculation based on the complex-contour integration 
technique. 

Consider the /i=0 case, i.e., the two lower degenerate bands are fully filled while the 
two upper degenerate bands are empty. Then a further A^-integral of Ci(k z ) gives the Hall 
conductivity a xy = —j- for d>0 and a xy = ^ — l) for d<0. The asymmetry of a xy 
between d>0 and d<0 will be explained below together with the behavior of the OM. When 
the local spins {Si} are inverted (which means that the spin chirality is also inverted), then 
the Hall conductivity changes its sign. 



IV. ORBITAL MAGNETIZATION AND ITS EFFECTS 

Now we turn to study the OM and its various effects. Without loss of generality, the 
present attention is only on the ^-component of the OM which is connected with Hall 
conductivity a xy . First, after a straightforward derivation, one obtains the fc-space orbital 
magnetic moment as follows 

« k > = 2^0" (26) 
See Eq. (|T9|) for F xy . One can see that the orbital magnetic moment is identical for upper 
and lower bands, while the Berry curvatures [Eq. ffl8]) ] for upper and lower bands differ by 
a sign. A comparison between Eq. fl26l) and Eq. ({TBI) gives an interesting relation for the 
present model 

m*(k) = ^(k)e nk . (27) 
11 




FIG. 4: (a) Zero-temperature orbital magnetization A4 Z and (b) Hall conductivity function 
of distortion d. The Fermi energy in (a) is chosen to be /zo=0.1, while in (b) the Fermi energy is 
chosen to be //o=0-0 which corresponds to case that only the lower bands are fully occupied. 

Given the expressions for m*(k) and Q^(k), the OM Ai z can now be systematically studied. 
At zero temperature, in particular, by substituting Eq. (127]) into Eq. PJ one finds an 
important relation 

M z = - h ^Y.\ ^n(k) (28) 

n 

_/fo 

e 

which indicates that the OM is proportional to the Hall conductivity with coefficient 
(— /io/e). Equation ( 128]) also holds at low temperature. Although this remarkable relation 
between the OM and the Hall conductivity is specific to the present Hamiltonian model, it 
definitely tells one that the topological ingredient in the OM may be faithfully mapped out 
through the Hall conductivity. If the band is partially filled, then after integrating by parts 
one finds that Eq. (1281) can be written as a pure Fermi-surface integral. In the integer filling 
case, on the other hand, the OM and Hall conductivity display a Fermi-sea feature. Figure 4 
shows M. z and function of the distortion d. The d>0 and d<0 cases are asymmetric 

by the observation that a xy is quantized for d>0 while non-quantized for d<0. To under- 
stand this asymmetry one may treat the distortion d as a control parameter of the 3D band 
structure and start from d=0, at which the upper and lower bands are degenerate along 
the BZ edge, i.e., ( y k x =±n,k y =±ii,k z ), (k x =,k y =±ii,k z =±7i), (k x =±n, k y , k z =±n) . When 
d is varied from d=0 to d>0, then these initial degenerate points completely split into two 
groups of "Dirac point" singularities, which play the role of positive and negative monopole 
sources, respectively [see, as an example, k and k' points in Fig. 2]. The upper and lower 
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bands are now tightly coupled by a series of "Berry flux loops". Along each loop Berry 
curvature flux 2tt passes from the lower bands to the upper bands through one Dirac point, 
then returns through the other corresponding one. The positive (negative) monopoles of the 
lower bands and the negative (positive) monopoles of the upper bands may recombine by a 
relative displacement of a primitive reciprocal lattice vector G. In the present cubic lattice 
the G is along one selective fc a -axis (a = x, y, z) with amplitude 2ir (the lattice constant a 
has been scaled to be unity). This means that during the "Dirac point" splitting process, 
the individual Chern invariant (the 3D generalization of 2D Chern number) for the lower 
and upper bands changes by =fG, respectively, while their sum is conserved. As a result, 
the Hall conductivity a xy for lower/upper bands is quantized to be a product of \ and the 
following Chern invariant 

i / d 3 kW n (k) = -^C n G z . (29) 

When d is varied from d=0 to d<0, equation (1291) breaks down, since the 2D Chern number 
C n is now fc-dependent by the presence of the additional two (3+l)D Weyl fermions at 
(k x , k y , k z ) = ±D(1, 1, 1). In this case, although the gap is opened at Dirac points along the 
BZ edge, the emergence of new Weyl fermions contributes non-zero density of states in the 
gap. It is then straightforward to repeat the integration in Eq. ( 1291) by parts to expose the 
non-quantized part of the 3D intrinsic Hall conductivity as a Fermi surface property. 

In the half filled case, i.e., when the chemical potential is in the Mott gap, the Hall 
conductivity vanishes due to the cancellation of the Chern invariants of the upper and lower 
bands. As a consequence, the OM in Eq. ( l29l) also vanishes. This result is prominently 



different from that in Ref. 



3, 



in which a finite OM with amplitude smoothly varied 



with distortion d was given via a tight-binding calculation. Ref. [16( employed the Kubo- 
Streda formula in calculating the OM. As we have shown above, however, the Kubo-Streda 
formula and the semiclassical formula give the same expression for the intrinsic OM. So no 
discrepancy is expected to occur between the two approaches. The non-zero OM in Ref. [3] 
in the Mott gap was ascribed by the authors to be a bulk Fermi sea property. However, the 
present explicit relation Eq. ( 1291) together with Haldane's argument on the metallic AHE 
shows that the non-quantized part of the OM in the present model is a Fermi-surface 
Berry phase effect, while the quantized part is completely determined by the topology of 
the filling bands. 
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FIG. 5: (Color online), (a) Orbital magnetization A4 Z , and (b) its two components M' (red curve) 
and Mq (blue curve) as a function of the electron chemical potential n for distortion d=0.2. The 
shaded area is the gap between the lower and upper bands. To suppress the divergence at band/gap 
contact, we have used the temperature of /cbT=0.05. 

Figure 5(a) shows the M. z as a function of the electron chemical potential /i for the 
distortion d=0.2. One can see that initially the OM rapidly decreases as the filling of the 
lower bands increases, arriving at a minimum at //=— 0.4, a value corresponding to the top 
of the lower band. Then, as the chemical potential continues to vary in the gap [shaded 
region in Fig. 5(a)] between the lower and upper bands, the OM goes up and increases as 
a linear function of [i. This linear relationship in the insulating region is explicit from Eq. 
( |28|) . When the chemical potential touches the bottom of the upper band, then the linear 
increase in A4 Z suddenly stops and the OM rapidly decreases again by going the chemical 
potential through the upper bands. The turning behavior at the band/gap contacts becomes 
numerically divergent at fcsT=0. This discontinuity is due to the singular behavior of fi n (k) 
at the BZ edge point k=ko, which will play its role when the fc-integral is over the entire 
BZ. 

The distinct behavior of the OM in the metallic and insulating regions, as shown in 
Fig. 5(a), reflects the different roles that its two components M* and Mq play hi these 
two regions. To see this, we show in Fig. 5(b) M* (blue curve) and (red curve) as a 
function of [i. One can see that and M^j oppose each other, which implies that they 
are carried by opposite-circulating currents. Also one can see that in the band insulating 
regime, the conventional term M* keeps a constant during variation of fi. This behavior is 
due to the fact that the upper limit of the k- integral of m n (k) is invariant as the chemical 
potential varies in the band gap. In the metallic region, however, since the occupied states 
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varies with the chemical potential //, thus M c also varies with fi, resulting in a decreasing 
slope shown in Fig. 5(b). The Berry phase term also displays different features in the 
insulating and metallic regions. In the insulating region, Mjj linearly increases with /i, as 
is expected from Eq. (|28[) . In the metallic region, however, this term keeps a constant with 
the amplitude sensitively depending on the topological property of the band in which the 
chemical potential is located. On the whole it reveals in Fig. 3 that the metallic behavior of 
the OM is dominated by its conventional term M c , while in the insulating regime the Berry 
phase term Mq comes to play a main role in determining the behavior of the OM. 

The above separate discussion of M c and can be transferred to study the anomalous 
Nernst effect (ANE). The relation between the OM and ANE has been recently found 
To discuss the transport measurement, it is important to discount the contribution from the 
magnetization current, a point which has attracted much discussion in the past. Cooper et 
al. [3l| have argued that the magnetization current cannot be measured by conventional 
transport experiments. Xiao et al. Q| have adopted this point and built up a remarkable 
picture that the conventional orbital magnetic moment M c does not contribute to the trans- 
port current, while the Berry phase term in Eq. ([31) directly enters and therefore modifies 
the intrinsic Hall current as follows 

p 2 r d 3 k 

j«= - ^ E >< E J J^M^ k)0>(k)-V x Mn(r), (30) 

In the case of uniform temperature and chemical potential, obviously, the second term is 
zero and the Hall effect of the distorted fee lattice is featured by the first term in Eq. (I30I) . 
In the following, however, we turn to study another situation, where the driving force is not 
provided by the electric field. Instead, it is provided by a statistical force, i.e., the gradient of 
temperature T. In this case, Eqs. (130]) and ([3]) give the expression of intrinsic thermoelectric 
Hall current as j x = ot xy {— V y T), where the anomalous Nernst conductivity a xy is given by 

1 e \ -s f d 3 k . . 

^y = TH^jj2^ nn (31) 



x 



[(enk - H) fn + k B T In (l + e^"^)] . 



Figure 6 shows a xy of the distorted fee lattice as a function of the chemical potential for 
d = 0.2 and /c#T=0.05. One can see that the ANE disappears in the insulating regions, 
and when scanning fi through the contacts between the bands and gaps, there will appear 
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FIG. 6: The intrinsic anomalous Nernst conductivity a xy as a function of the electron chemical 
potential /i for d=Q.2 and fesT=0.05. The shaded area is the gap between the lower and upper 
bands. 

peaks and valleys. Remarkably, a similar peak-valley structure was also found by the recent 
first-principles calculations in CuCr2Se4 compound [3]. The ANE of this compound was 
recently measured by Lee et al. [6j as a function of Br doping which substitutes Se in the 
compound and changes the chemical potential \i. Due to the scarce data available, until 



now the peak- valley structure of a xy revealed in Fig. 6 and in Ref. js] has not been found 
in experiment, and more direct experimental results are needed for quantitative comparison 
with the theoretical results. Interestingly, the expression for a xy can be simplified at low 
temperature as the Mott relation 



7T 



2 k%T da xy {/i ) 



^ = -t-t~^- (32) 

Thus one can see that the low-temperature nonzero ANE is a Fermi-surface Berry phase 
effect. Another unique feature of a xy is its linear dependence of temperature. 



V. CONCLUSION 



In summary, after pointing out the equivalence of the semiclassical approach and the 
quantum Kubo-Streda formula in description of the intrinsic orbital magnetization, we have 
theoretically studied the properties of the OM in antiferromagnets on the distorted 3D fee 
lattice. The distortion parameter d in the fee lattice produces nonzero 2D Chern number 
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and results in profound effects on the OM properties. An explicit relation between the 
OM and the Hall conductivity in this system has been derived. According to this relation 
we have found that the OM vanishes when the electron chemical potential is lies in the 



Mott gap, which is in contrast with the results in Ref. 16j. We have shown that the two 
parts M c and Mjj in the OM oppose each other, and yield the paramagnetic and diamagnetic 
responses, respectively. In particular, due to its Fermi-sea topological property, the magnetic 
susceptibility of remains to be a nonzero constant when the Fermi energy is located in 
the energy gap. It has been further shown that the OM displays distinct behavior in the 
metallic and Chern-insulating regions, because of different roles M c and Mo play in these 
two regions. The anomalous Nernst conductivity has been studied, which displays a peak- 
valley structure as a function of the electron chemical potential. We expect that these results 
will be experimentally verified in the spin-frustrated systems. 
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